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We discuss the extraction of negative-parity baryon masses from lattice QCD calculations. The mass of the 
lowest-lying negative-parity J = 1/2“ state is computed in quenched lattice QCD using an Cl(a)-improved clover 
fermion action, and a splitting found with the nucleon mass. The calculation is performed on two lattice volumes, 
and three lattice spacings enabling a study of both hnite-volume and finite-lattice-spacing uncertainties. A 
measurement of the hrst excited radial excitation of the nucleon finds a mass considerably larger than that of the 
negative-parity ground state, in accord with other lattice determinations but in disagreement with experiment. 
Results are also presented for the lightest negative-parity 7 = 3/2 state. 


1. INTRODUCTION 

The N* spectrum exhibits many features that 
are emblematic of QCD; that they contain three 
valence quarks is a manifestation of SU(3) sym¬ 
metry, while the linearly rising Regge trajectories 
are evidence for a ffux-tube-like confining force. 
Finally, the fine and hyperfine structure reveals 
details of the interquark interactions. Thus the 
study of the spectrum of excited nucleon reso¬ 
nances has always been a vital component of the 
experimental hadronic physics program, as typi¬ 
fied by the Hall B experimental effort at Jefferson 
Laboratory. 

The motivation for studying the N* spectrum 


in many ways mirrors that for the hybrids, which 
also provide an arena in which to explore flux 
tubes. One distinction is that, in the case of nu¬ 
cleons, all spin-parity combinations are accessible 
to the quark model, and thus there are no nu¬ 
cleon “spin exotics”. However, both hybrids and 
excited nucleons share the feature that they are 
sensitive to the presence of excited glue. 

The observed excited nucleon spectrum al¬ 
ready poses many questions. There are “miss¬ 
ing resonances”, predicted by the quark model 
but not yet observed experimentally. Most strik¬ 
ingly, there is the so-called Roper resonance, and 
the lightest negative-parity A state, the A(1405) 
with anomalously light masses suggesting that 
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they may be “molecular” states, rather than true 
three-quark resonances. Lattice gauge theory cal¬ 
culations have a vital role to play in resolving 
these questions. 

In an attempt to address these issues, there 
have been several recent calculations of the 
lowest-lying excited nucleon masses, emphasising 
in particular the extent to which the parity part¬ 
ners are accessible to lattice calculation. The 
first calculation employed the highly-improved 
D 234 fermion action|j^,||, whilst a second calcu¬ 
lation employed domain-wall fermions [^|J^]. Both 
calculations exhibited a clear splitting between 
the masses of the A^^/^+(937) and A^^/^“(1535) 
states, with an indication that the lowest-lying 
radial excitation of the nucleon has a mass consid¬ 
erably larger than that of the A^*(1535), casting 
doubt on the nature of the Roper resonance. 

In this talk, we present a calculation of the low¬ 
est lying negative parity nucleon mass using an 
O(a)-improved Sheikholeslami-Wohlert (SW), or 
clover, fermion action; preliminary results were 
presented in ref. 1^, and a comprehensive analy¬ 
sis appeared in ref. . By choosing the coefficient 
of the improvement term appropriately, all 0 {a) 
discretisation uncertainties are removed, ensur¬ 
ing that the continuum limit is approached with 
a rate proportional to a^. For a subset of our lat¬ 
tices, masses of the first radial excitation of the 
nucleon and of the parity partner of the A are 
presented. 

The calculation of the excited nucleon spec¬ 
trum places particularly heavy demands on lat¬ 
tice spectroscopy. The excited nucleon states are 
expected to be large; the size of a state is ex¬ 
pected to double with each increase in orbital an¬ 
gular momentum. Thus a lattice study of the ex¬ 
cited nucleon spectrum requires large lattice vol¬ 
umes, with correspondingly large computational 
requirements. Furthermore, the states are rela¬ 
tively massive, requiring a fine lattice spacing, at 
least in the temporal direction. These require¬ 
ments could be satisfied with much greater econ¬ 
omy using the clover fermion action than using 
the domain-wall or overlap formulation. Thus it 
is important to establish that the negative par¬ 
ity states are indeed accessible to calculations us¬ 
ing the clover action. Finally, by comparing the 


masses obtained using the clover action with a 
calculation, at a single quark mass, using the Wil¬ 
son fermion action, we also gain insight into the 
nature of the interaction responsible for the split¬ 
ting in the parity doublet. 

The layout of the remainder of this talk is as fol¬ 
lows. In the next section, I shall describe the con¬ 
struction of hadronic operators, and in particular 
the projection onto positive- and negative-parity 
states, and summarises other computational de¬ 
tails. Section H contains the results for the masses 
of the lowest-lying positive- and negative-parity 
nucleon using the SW fermion action, and for 
a subset of the lattices the masses obtained for 
the Roper resonance. The talk concludes with a 
discussion, a comparison with other lattice and 
phenomenological calculations, and prospects for 
future studies. 

2. SIMULATION DETAILS 

2.1. Baryon operators 

Historically, there has been relatively lit¬ 
tle study of the negative-parity baryon states, 
though there had been numerous computations 
of the nucleon radial excitations. It is therefore 
useful to review the role that parity plays in the 
construction of baryon operators, and how this 
role is manifest in the case of lattice calculations. 
We illustrate the discussion through considera¬ 
tion of the usual nucleon interpolation operators: 



= e^jk{ufC-f5dj)uk, 

( 1 ) 



( 2 ) 


= e^jk{ufC"f4j5dj)uk. 

( 3 ) 


These operators have an overlap with particle 
states of both positive and negative parities. In 
order to construct operators of definite parity, we 
introduce the parity-projection operator: 

P=(I±74), (4) 

and construct the correlators 

C±it) = ^(0|TiV(f,t)(l ±74)A(0)|0 >, (5) 

X 

where iV is a baryon interpolating operator. On 
a periodic or anti-periodic lattice, the two time 
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orderings result in forward propagating states of 
positive (negative) parity, and backward propa¬ 
gating states which are anti-particles of negative 
(positive) parity states if the projection opera¬ 
tor is chosen with positive (negative) sign. Thus 
the best delineation that may be constructed is 
between a forward propagating particle, and the 
backward propagating anti-particle of the parity 
partner. At large distances, where t ^ 1 and 
Nt — t ^ 1 , the correlators assume the time be¬ 
haviour 

C+it) ( 6 ) 

C-it) -f (7) 

where mf and m~ are the lightest positive- and 
negative-parity masses respectively. 

For the case of the nucleon ground state, the 
“diquark” part of both Ni and couples upper 
(large) spinor components, while that in N 2 in¬ 
volves both an upper and lower spinor component 
and thus vanishes in the non-relativistic limit [Q. 
Thus our expectation is that the operators Ni, 
N 3 should give a better overlap with the positive- 
parity ground state nucleon compared with oper¬ 
ator N 2 , and this is confirmed in lattice simula¬ 
tions. For some of the lattices used in the calcula¬ 
tion, the positive-parity nucleon mass is obtained 
using the non-relativistic quark operators Q| , de¬ 
fined by 

i’^ + = ip^{i + 74)- ( 8 ) 

corresponding to the “large” components in the 
non-relativistic limit; this results in a factor of 
two reduction in computational effort at the ex¬ 
pense of some loss of statistics. 

2.2. Simulation Details 

The calculation is performed in the quenched 
approximation to QCD. The standard Wilson 
gluon action is employed, whilst the quark prop¬ 
agators are computed using the Sheikholeslami- 
Wohlert clover fermion action, using the non- 
perturbative coefficient of the clover term deter¬ 
mined in refs. Q and j|]. Thus the hadron masses 
are free of all 0 {a^) discretisation errors. 

The propagators are computed using both lo¬ 
cal sources, and spatially extended “smeared” 


sources, using either the fuzzed-source method 
(F) or Jacobi-smeared sources (J) as described 
in references 0 and [|ll[ respectively. The sim¬ 
ulation parameters are summarized in Tables |^. 
The errors on the fitted masses are computed us¬ 
ing a bootstrap procedure. Different numbers of 
configurations are used at different quark masses 
at some of our /3 values, precluding the use of cor¬ 
related hts in some of the chiral extrapolations. 
Thus a simple uncorrelated fit is performed, 
with the uncertainties computed from the varia¬ 
tion in the x^- Calculations of the light hadron 
spectrum using these lattices were presented in 
references 0 and 00 - 

The clover term removes the leading chiral- 
symmetry-breaking effects from hadron masses at 
finite a. Since the lack of degeneracy between the 
positive- and negative-parity baryon states is a 
consequence of the spontaneous breaking of chi¬ 
ral symmetry, we also compare the measurements 
of the mass splitting at a single quark mass with 
those obtained using the Wilson fermion action. 

3. RESULTS 

The masses of the lowest-lying and 

states are obtained from a simultaneous, 
four-parameter fit to the positive- and negative- 
parity correlators C+{t) and C-{t) constructed 
from the “good” nucleon operator Ni of eqn. |^. A 
clear signal for the mass of the lightest negative- 
parity state notably higher than that for the 
positive-parity state is observed, and is illustrated 
for one of our ensembles in Figure ^ 

For the chiral extrapolation of the hadron 
masses, we adopt the ansatz 

{amxf = {aMxf + 62 ( 0771 ^)^ (9) 

where we use upper-case letters to denote masses 
obtained in the chiral limit and X is either 
or . We include data at different volumes, 

but at the same in the chiral extrapolations, 
but treat the fuzzed and Jacobi-smeared data in¬ 
dependently; they give slightly different masses, 
reflecting the different contributions of the ex¬ 
cited states in the two cases. Note that the fit 
to {amxY' is formally the same as one to amx 
at small pion mass, but in practice provides a no- 
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Table 1 

The parameters of the lattices used in the calculation. The labels J and F refer to use of Jacobi and 
“fuzzed” quark sources respectively. 


p 

^SW 

-T 

L [fm] 

K 

Smearing 

6.4 

1.57 

32^ • 48 

1.6 

0.1313, 0.1323, 0,1330, 0.1338, 0.1346, 0.1350 

J 

6.2 

1.61 

243 • 48 

1.6 

0.1346, 0.1351, 0.1353 

F 



243 • 48 

1.6 

0.1333, 0.1339, 0.1344, 0.1349, 0.1352 

J 



323.04 

2.1 

0.1352, 0.1353, 0.13555 

J 

6.0 

1.76 

163 • 48 

1.5 

0.13344, 0.13417, 0.13455 

F 



163 • 32 

1.5 

0.1324, 0.1333, 0.1338, 0.1342 

J 



243 • 32 

2.2 

0.1342, 0.1346, 0.1348 

J 


Lattice sizes in physical units are quoted using ro to set the scale |l^. 



Figure 1. The effective masses for the 
channel (circles) and the channel (dia¬ 

monds) at /3 = 6.2, with k = 0.1351, using the 
interpolating operator of eqn ^ Also shown 
as the bursts is the effective mass of the positive- 
parity state determined using the interpolating 
operator A^ 2 - 


tably better description of the data. 

It is known that there are non-analytic terms 
contributing to the chiral extrapolation of physi¬ 
cal quantities, and the contribution of the non- 
analytic terms arising from the pion-induced 
baryon self energies has been explored in refer¬ 
ence [ p^ . In the quenched approximation, the 
leading, non-analytic terms can assume a differ¬ 
ent form. From quenched chiral perturbation the¬ 
ory, one expects leading non-analytic terms which 
are linear in |^; the coefficient of the lead¬ 
ing term is predicted to be negative. We therefore 
also attempted to fit the data at /3 = 6.4 to the 
form 

{amx)^ = {aMxY + bi(amTr) + 62 ( 0771 ^)^. ( 10 ) 

Fitting the state we found a positive value 

of 61 , though with a very large error that would 
still accommodate a negative value. Thus it is 
unclear whether eqn. ^ provides a reliable form 
with which to extrapolate data obtained with 
quark masses around that of the strange to the 
chiral limit, emphasising once again the need for 
data at the smaller quark masses at which the 
pion cloud emerges. The chirally extrapolated 
masses obtained using this procedure, with bi 
unconstrained, differ from those obtained using 
eqn. ^ by around 5%. In view of the difficulties 
discussed above, and because on some of our en¬ 
sembles we have insufficient data points to per¬ 
form a three-parameter fit, we quote our final re¬ 
sults from fits to eqn. An example of such a 
chiral extrapolation for the data at /3 = 6.4 is 





















Figure 2. The masses in lattice units of the 
lowest-lying positive- and negative-parity nucle¬ 
ons at /3 = 6.4. The curves are from fits to 
using eqn. 


shown in Figure 

Ideally, the forms eqn. ^ and should be ap¬ 
plied to the baryons masses obtained in the in¬ 
finite volume limit. Indeed, the talk by Stuart 
Wright at this workshop emphasised the impor¬ 
tance of correctly accounting for finite-volume ef¬ 
fects in chiral extrapolations |Q. Our data do 
not admit this procedure, but an investigation of 
the masses of the negative-parity state obtained 
on the smaller and larger of our lattices at /3 = 6.0 
and 6.2 suggests that the masses of the negative- 
parity state on the larger lattices are smaller than 
those on the smaller lattices by around 5%. 

In order to look at the discretisation uncertain¬ 
ties in our data, we show in Figure § the masses 
in units of tq against the where rg = 0.5 fm 

is the hadronic scale |l^Jl4|| . 

On a subset of our lattices, we have extracted 
the ground-state positive- and negative-parity 
masses using the operator N 2 of eqn. ||. The ef¬ 
fective mass in the positive-parity channel using 
this operator is shown in Figure |^; it indicates an 
effective ground state mass considerably higher 
than that of the lightest negative-parity state. 
The signal for both the positive- and negative- 





Figure 3. The masses of the lowest-lying positive- 
and negative-parity baryons in units of 
against in units of r^. The lines are linear 
fits in a^/rg to the positive- and negative-parity 
baryon masses. Also shown are the physical val¬ 
ues. 


parity ground states degrades appreciably with 
decreasing pion mass, as can be seen in Figure ^ 
where we show the extracted lowest-lying masses 
of both parities using each of the A^i and N 2 in¬ 
terpolating operators. A better determination of 
the first radial excitation would require the mea¬ 
surement of the cross-correlator between 7Vi and 
A^ 2 , as done in reference |^]. However, the conclu¬ 
sion that the operator N 2 has a negligible overlap 
with the ground state nucleon at these values of 
the pseudoscalar mass seems inescapable. 

The measurement of heavier excitations is in¬ 
evitably noisier because of the reduced signal-to- 
noise ratio. On some of our lattices we have mea¬ 
sured the correlators for the I = 3/2 A, using an 
interpolating operator 

^ 3 / 2 . 1/2 ^ eijkiuJCj^Uj)uk. ( 11 ) 

This has an overlap onto both spin-3/2 and spin- 
1/2 states, but these can be distinguished using 
a suitable spin projection, and the two parities 
delineated as described above. The masses of the 
ground state of both positive and negative parity 
are shown in Figure ||. Further success at deter¬ 
mining the higher spin states has been made in a 
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Figure 4. The masses in lattice units of the 
lowest-lying (circles) and the (di¬ 

amonds) obtained using Ni, and of the lowest- 
lying (down triangle) and (up tri¬ 

angle) using ‘N 2 , obtained on the 24^ x 48 lattices 
at /3 = 6.2. 


recent determination of the masses ||2^ . 

4. CONCLUSIONS 

We have seen that the use of the clover- 
improved fermion action enables the delineation 
of the positive- and negative-parity excited state 
masses, in accord with calculations using highly 
improved Wilson, and domain-wall fermion ac¬ 
tions. On a subset of our lattices, we also com¬ 
puted propagators using the unimproved Wilson 
fermion action at a single quark mass, and found 
a parity mass splitting consistent with that of 
the clover fermion action. On reflection, this is 
hardly surprising; under SU{6) spin-flavour sym¬ 
metry, the low-lying positive-parity baryons can 
be assigned to an / = 0 multiplet, whilst the low- 
lying negative-parity baryons can be assigned to 
an 1 = 1 multiplet. Thus the mass-splitting can 
be likened to the P — S splitting in the meson 
sector, which is well-known to be faithfully re¬ 
produced using the Wilson fermion action. 

The striking feature common to all the lattice 


Figure 5. The masses in lattice units of the A (cir¬ 
cles) and its parity partner (diamonds) against 
the pseudoscalar mass obtained on the 24^ x 48 
lattices at /3 = 6.2 


calculations is the apparent inversion of the order¬ 
ing of mass of the lowest radial excitation of the 
nucleon, and the mass of the parity partner of the 
nucleon. All the calculations, using different glu- 
onic and fermionic actions, suggest a splitting be¬ 
tween the masses of the first positive-parity radial 
excitation and the lightest negative-parity state 
comparable to that between the masses of the 
negative-parity state and the nucleon. In con¬ 
trast, the observed lightest radial excitation of 
the nucleon is the Roper resonance with a mass 
of 1440MeV, whilst the lightest negative-parity 
state is the heavier A* (1535). We can see that 
lattice calculations of the spectrum are already 
posing interesting questions about the nature of 
the Roper resonance. 

I have not discussed the issue of the nature 
of the A(1405) in this talk. In the case of the 
positive-parity states, lattice calculations, such as 
that of ref. [^, suggest that the quark mass be¬ 
haviour of the baryon masses, in the range over 
which they are measured, is well described by 
terms corresponding to the “centres of gravity” of 
the quark masses within the hadron. This expec¬ 
tation has been confirmed in ref. 0 , where the 
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lightest negative-parity A resonance has a mass 
considerably in excess of 1405 MeV. It should 
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be noted, however, that a recent large Ng cal- 


culation has been able to accommodate a light 
A(1405) m. 

We have seen that lattice calculations are al¬ 
ready addressing some of the salient issues in 
baryon spectroscopy. The calculation of the ra¬ 
dial excitations and the higher spin states will 
present further challenges; the need to overcome 
the reduced signal-to-noise ratio afflicting heav¬ 
ier excitations, the construction of suitable lat¬ 
tice operators to isolate higher excitations, and 
the exploitation of improved data analysis meth¬ 
ods. As these issues are tackled, we will be able to 
tackle the more ambitious projects of transition 
form factors and hadronic decays. 
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